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The programs in this report are for use 
within the Department of the Navy, and they are 
presented without representation or warranty of 



any kind. 



A PROCEDURE FOR ESTIMATING AN OBJECT'S POSITION 



BASED ON TWO OR MORE BEARINGS WITH A 
PROGRAM FOR A TI-59 CALCULATOR 

I . Introduction 

A procedure for estimating an object's position with 
bearings taken on or from two or more stations is developed 
in Section IV of this report. In the development of the pro- 
cedure, the following things are assumed: The object and the 

stations are fixed on the surface of a flat earth and the 
position of each station is known. The error in the bearing 
taken on or from a station is a normal random variable with 
a known standard deviation e and a mean of zero (if bias 
exists, it is known and removed); and station bearing errors 
are independent. The user instructions for a TI-59 program 
to implement the procedure are given in Section II, and the 
program listing is given in Section III. 

As an example to illustrate a use of the program, suppose 
bearings are taken on an object from three stations (1, 2 and 3) 
as illustrated in Figure 1. Also, suppose that the assumptions 
stated above are satisfied and that an initial estimate of the 
object's position is made and that it is relatively near the 
object. This assumption is discussed in Section IV. 



1 



Station 3 




FIGURE 1. Geometry for the Example 



Let the measured bearings and bearing errors (standard 
deviations) be: 



= 35° 
e 2 = 35i° 
e 3 = 131° 

And let the ranges and bearings 




of the initial estimate be: 



2 



10,000 meters, 



38° 



r i - 



B i - 



= 15,000 meters, 



r 3 = 12,000 meters, 



e 2 - 346 ° 



6, = 127° . 



Use of the position estimation program with this data 
gives a final position estimate (fix) determined by: 



x = -512 meters 
y = - 75 meters 

where x is its East-West distance and y is its North-South 
distance from the initial position estimate. The East-West, 
North-South xy-coordinate system with its origin at the initial 
estimate is shown in Figure 2. So the final position estimate 
is 512 meters to the West and 75 meters to the South of the 
initial position estimate. 

North 

y 





Initial Estimate 


* Final Estimate 





x East 



FIGURE 2. The Location of the Final Position Estimate with 
Respect to the Initial Position Estimate. 
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Minimum area elliptical confidence regions for an object's 
position can also be found by using the TI-59 program. The 
centers of the regions are at the fix, and their axes lie along 
the x' and y' axes of the coordinate system obtained by 
rotating the East-West, North-South xy-coordinate system with 
origin at the fix through an angle Y . The angle Y is defined 
so that it is positive for a rotation in the counterclockwise 
direction. 

With the data from the above example, the program gives 
Y = -31°; so, the x' axis is directed 31° South of East. For 
a confidence region with minimum area and a confidence level of 
.9000, the ellipse bounding the region has a semi-major axis 
of 2064 meters, and a semi -minor axis of 1453 meters. The area 
of the region is 9.43 square killometers or 2.75 square nautical 
miles. The region is shown in Figure 3. 



North 




FIGURE 3. A .9000 Confidence Region for an Object's Position. 
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In the example discussed above, the position of the 
initial estimate is an input to the program. If this is 
not desirable, the program can be used to determine a po- 
sition for the initial estimate. The position is the in- 
tersection of the two bearing lines corresponding to the 
first two bearings entered in the program. Both options 
are illustrated in Section II. 

Since, in general, the smaller the bearing errors, 
the more likely that the initial estimate will be relatively 
near the object; small bearing errors can be considered to 
be a condition on the use of the procedure. 

Note, if the length of the base line joining the 
first two stations is small enough and their bearing errors 
are large enough, observed bearing lines from the two 
stations may not intersect. If they do not intersect, the 
initial estimate determined by the program will be at the 
intersection of the reciprocal bearing lines, and a gross 
error can result. 
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II. User Instructions 



The TI-59 program to which the user instructions in this 
section apply can be used to calculate the quantities described 
in Section I. 

The program requires the following inputs: 

1. the observed bearing from or on an object for two or more 
stations ; 

2. station positions relative to a reference position; and 

3. the bearing error (standard deviation) for each observed 
bearing . 

Station positions can be specified in either of two ways. 

In the first way. Mode A, each station's position is specified 
in terms of its bearing a and its range p from a reference 
position. In the second way. Mode B, each station's position 
is specified in terms of its East-West distance x (plus for 
East) and its North-South distance (plus for North) from a 
reference position. The reference position can be any convenient 
location. For example, if it were at a station, then for 
that station a = 0 and p = 0 or x = 0 and y = 0. 

The program also requires an initial estimate of the object's 
position. The user has two options: 

1. Let the program provide an estimate, or 

2. Provide one with the input data. 

For Option l,the initial estimate is at the intersection of the 
bearing lines determined by the first two observed bearings entered 
into the program. For this reason, if this option is chosen, the 
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first and second groups of data entered should correspond to 
the two stations estimated to have the smallest products r^e^. 
Although in this option the reference position cannot be at 
the initial estimate, it can be at one of the stations. If 
only two stations are involved, the final estimate is at the 
intersections of the bearing lines. (If the second option of 
either mode is used with an initial estimate which is not at 
the intersection of the two bearing lines, the coordinates of 
the final estimate will differ from coordinates of the inter- 
section to the degree of the approximations involved in the 
estimation procedure.) 

Two ways of providing confidence (probability) region data 
are available. In the first way, Mode C, the confidence 
(probability) p is specified. In the second way. Mode D, the 
multiplier k is specified where ka~ , and ka~ , are the 
semi-axes of the bounding ellipse. 

The values of various quantities calculated by the program 
are either stored in registers or appear in the display. If a 
PC-100A printer is used, some of these values will be printed. 

The location of calculated values and the printing format is 
given after the user instructions. Those quantities which are 
not described in Section I are described below in the User In- 
structions or in Section IV. 

All angles required or calculated by the program are 
in decimal degrees. 
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Step 



Instructions 



Enter Press Display 



1. If the calculator has been in use and 
flags have been set or the memory 
repartitioned, turn the calculator off 
and then on. 

^ 2. Read Side 1 and Side 2 of Card 1. 

3 . Read Side 3 of Card 2. 



MODE A: Station Locations Specified in 

Terms of Bearing and Range from 
a Reference Point. 

4a. If the initial position estimate will 
be determined by the program, go to 





Step 


7a. 


See the 


note on Page 10. 






5a. 


Enter 


the 


initial 


estimated bearing. 


a* 


A 1 


6a. 


Enter 


the 


initial 


estimate's range. 


P* 


R/S 


7a. 


Enter 


the 


measured bearing on the 


0 i 


A 



object from a station or the reciprocal 
of the measured bearing on a station 
from the object. 



8a . 


Enter 


the 


station's bearing. 


a i 


R/S 




9a. 


Enter 


the 


station's range. 


Pi 


R/S 




10a. 


Enter 


the 


bearing error. 


e . 

l 


R/S 


i 



11a. Repeat Steps 7a, 8a, 9a and 10a for all 
stations. The number of repetitions i 
appears in the display after Step 10a. 
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Step 



Instructions 



Enter Press 



Display 



^ MODE B: Station Locations Specified in Terms 

of East-West Distance and North-South 
Distance from a Reference Point. 



v 4b. If the initial position estimate will 
be determined by the program, go to 
Step 7b. See the note on Page 10. 



5b. 


Enter the initial estimate's East-West 
distance . 


X* 


B' 




6b. 


Enter the initial estimate's North- 
South distance. 


y* 


R/S 




v 7b. 


Enter the measured bearing on the 
object from a station or the 
reciprocal of the measured bearing 
on a station from the object. 


e i 


B 




*8b. 


Enter the station's East-West 
distance. 


x i 


R/S 




^9b. 


Enter the station's North-South 
distance . 


Vi 


R/S 




1.10b. 


Enter the bearing error. 


e i 


R/S 


i 



* lib. Repeat Steps 7b, 8b, 9b and 10b for all 
stations. The number of repetitions i 
appears in the display after Step 10b. 



BOTH MODES 

*12. Calculate the East-West distance, R/S 

the North-South distance, the 
bearing and the range of the 
position estimate relative to the 
reference position. Also calculate 
the rotation angle y , and the 

standard deviations O^, and 0^,. 

x 1 y 1 

To include additional bearing measure- 
ments after this calculation, go to 
Step 18. 
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Step 



Instructions 



Enter 



Press 



Display 



13. 


For confidence (probability) region 
calculations, go to Step 14 if the 
confidence (probability) for the 
region is specified. If k is 
specified where ko^ ( and ko Ai 

X y ' 

are the semi-axes of the bounding ellipse 
with the larger the major axis, go 
to Step 16. 


14. 


Enter p, the confidence level p C Area 

(probability) and calculate k, 

ko A , , kn and the area of the region. 

X ' y ' 

(The area units correspond to the 
distance units used.) 


^15. 


For a different value of p, go to 
Step 14. 


16. 


Enter k and calculate the confidence k D Area 

level (probability) p, ko^ # , ko A1 and 

x 1 y 1 

the area of the region. (The area 

units correspond to the distance units 
used. ) 


17. 


For a different value of k , go to 
Step 16. 


18. 


To include an additional bearing measure- 
ment from either a new or old station, 
go to Step 7a if using Mode A or 
Step 7b if using Mode B. 



NOTE: 


If a data entry error occurs in either mode, 
press RST and then use the following pro- 
cedure: For Option 1, return to Step 7 

and repeat all data entries. For Option 2, 
return to Step 5 and repeat all data entries. 

Also, if a position estimate is to be de- 
termined for a new object position or if 
a new mode is to be used, follow this 
instruction . 
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NOTES : 



a) The program printing format is given below: 

For the initial data, i = 1, 2, ... , n with one space 
between groups : 



Mode A 
a* 

P* 




e . 

1 



initial estimate 
if provided 



Mode B 
x* 

y* 

e . 

i 

X . 

1 

y i 

e . 

i 



The format for the calculated position data is: 

x 

y 



a 

P 

Y 
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For the confidence (probability) region portion of the program 
the format is after pressing either C or D: 



P 

k 



x ' 


semi-axis 


ko A , 

y' 


semi-axis 


Area 





b) The following data is 


stored in the indicated registers 


Data 


Registers 


X* 


R38 


y* 


R39 


Y 


R29 


X 


R30 


y 


R31 


a 


R32 


P 


R33 


a - 

x ' 


R16 


a 

y 


R17 



P 


R14 


k 


R15 


ka~ , 


R18 


ka~ , 

y’ 


R19 

12 



R19 



Four data tapes for a sample problem are given below. 



Distance units have not been specified, but they could be 
meters for example. Angles are in degrees. Option 1 (initial 
estimate not provided) for Mode A and Mode B is indicated by 
A and by B and Option 2 (initial estimate provided) is indicated 
by A ' and B ' . 

For each mode and each option, the input data are 
indicated. The data determine the relative locations of three 
stations as well as the observed bearing of an object from 
each station. 

For A and B, the reference location is at Station 1 
and the initial position estimate (determined by the program) 
is at the intersection of the bearing lines for Station 1 and 
Station 2. 

The intersection has coordinates x* = 906.4853528 
and y* = 17296.77092 with respect to Station 1. 

For A' and B', both the initial estimate and the 
reference location are at the intersection of the bearing lines, 
so a* = 0 and p* = 0 and x* = 0 and y* = 0 . 

The data for A, B, A' and B' are all equivalent, 
and each solution gives the same data for a confidence (prob- 
ability) region calculation. A tape with confidence (prob- 
ability) region results for both Mode C and Mode D which 
correspond to A, B, A' and B' is given with the first four 
data tapes. 
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■J s' "= = '3 

1 6462 
1= 995 
19472 




A 



A’ 



0. 

U c 

4. 



e i 

“l 

Pi 

e. 




i 0000. 



a 2 

P2 

e 2 



303. e 3 

o t ot 2 

14000. P 3 



S78933 x 
. 71223 y 
471725 a 
.70157 p 



392245 Y 
963755 a ~, 

080777 a* 

y’ 



0. a* 

0. P* 

. .3- 01 

i 8 2 . ai 

17320.50808 p x 



33. e 2 
213. a 2 
20000. p 2 
3. e 2 

303. 0 3 

129.5867755 « 3 
8717. .-97886 P 3 

8. e 3 

- 3 ’ :• 2 . 8 97 45 6 7 x 

” 8 4 . U 5 8 6 8 8 5 v 
201.7584019 <* 

898. 03931 ii p 

-7. 325392259 y 
787. 3663757 a 
1233.080776 a 
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^<>X> 



B 



C or D 



B' 



U. 

0. 

4. 



-9986. 295348 

c- o o --i cqCi; • j ^ 



°2 

x 2 

*2 

e 2 



303. 
7624. 94649 
11741. 38795 



x 3 
y 3 

e 3 



573. 5878927 x 
1 646 * 2 . 71223 y 
1.995471723 a 
16472. 7 U 1 5 7 p 



-7. 325392245 
787. 3663756 
1233. 080777 



Y- 




y' 



0. x * 
0. y* 

3. 0JL 
•906. 4853528 x ± 
* 1 i 2 9 6 . 7 7 0 9 2 y ]_ 

4 • ei 



- i u 8 9 a . r* y u 7 X 2 
1 6 7 7 y . 4 i i y 6 y 



0 2 

2 

2 



303. 

671 y . 4 6 1 1 y 7 
*5555. 382969 



0 3 

x 3 

y 3 

e 3 



• 332 . 8974589 x 
■ 8 8 4 * U 5 8 6 9 U 6 * y 
201. 7584019 a 
898. 0393184 p 



i . •-« C • J O 7 cl ? 

787. 3663755 
t 2 *J 3 « U 8 LI r r !*’ 



Y 

a. 



y' 



0. 9 


P 


t 1 4 5; y h: r s f j 2 h 


k 


I 6 8 9 . t 149 3 


kcu 


'V J. £ 1 A Q .-j cr .J 

i- v -• c* 8 i. h- r 4 j •-+ 


, X 1 

kO ; ; , 


14046 564. 97 


y 

Area 


e O 6 ^ h h 4 7 J rj : ;-! 


P 


e 


k 


1574. 732751 


k ° 


2466. 161553 


ka§, 


1 2 2 U U 5 1 7 . 6 


Area 
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The following data tape illustrates the effects of 
using only the bearings for the first two stations. The tape 
is for Mode A, Option 1. The values stored in Registers 38 
and 39 (the x and y coordinates of the intersection of the 
bearing lines from Station 1 and Station 2 with reference, 
in this case, to Station 1) are also listed on the tape (as 
well as given above). And, as can be seen, the initial 
estimate and final estimate correspond. 

The data tape also illustrates the use of additional 
bearing data to revise a position estimate. The data for 
Station 3 printed after the first confidence (probability) 
region calculation results was entered by again repeating 
Step 11a, and the remaining results were obtained by next 
repeating Step 12 and then Steps 14 and 16. Note, these 
results are the same as the corresponding results for Mode A 
on page 14. 
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A 



0. 

0. 

4. 



6i 

a l 

Pi 

e. 



•7. 325392245 
787. 3663755 
1233. 080777 



Y 

a -, 

x 

a y» 



ti r o c 

1 0000 . 



906. 4853528 
17296. 77092 

•— 1 C 

17320. 50808 



5750198 
. 886822 
3092. 663848 



- 20 . 3 
8 1 8 



02 

a 2 

P2 



x 

y 

a 

P 



Y 

a 

A f 

a5, 

y 



0. 9 

2. 145966026 
1689. 661492 
c. fc« 4 b c i 4 y 4 5 u 
14046364. 98 



. o 6 4 6 6 4 7 i 6 y 

e 

-* cr “7 A — ~i ~7 ET H 

1 j i - 4 . i £ i j i 

2466 . 161 5 5 4 
12200517. 6 



P 

k 

ka -. 

ka* 

y 

Area 



P 

k 

kO/v, 

i x 

ka-. 

Area 



0. 9 

145966 Li 2 fc» 
1757. 303299 
6 6 3 b g 75154 9 
3 6 6 3 973 U .91 



P 

k 

ka„, 

ka^, 

Area 



. 8 64664716 8 


p 


*~l 
U* C 


k 


1 6 y 7 g 77 3 64 4 


ka -, 


6 1 y 5 « y 2 7 6 9 6 


ka£, 


1 y 2 4 y 5 1 *’ « 2 2 


Area 



906. 4 853528 R38 

17296.77092 R39 



303. 
o o . 
14000. 



ou 



573. 5878933 
16462. 71223 
1. 995471725 
16472. 70157 



X 

y 

a 

P 
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To obtain the results given in Section I, use A’ 
and take the reference position at the initial estimate 
(a* = 0, p* = 0). Then = 218°, a 2 = 166° and 

= 307°. The data tape for the calculation is given below. 
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A’ 



0. 

0. 


a* 

P* 


35. 
218. 
1 0000. 

4. 


01 

a i 

Pi 

e i 


351. 


02 


166. 
1 5000. 

( . 


a2 

P 2 

e 2 


1 3 1 . 

307. 

12000. 

cr 

o 


0 3 

a 3 

p 3 

e 3 


-511. 961856 
-75. 43753333 
261. 617789 
517. 4398687 


X 

y 

a 

P 


-31. 23492683 
677. 2632305 
96 i . 6888632 


Y 

a. , 

x' 

a ^ 

y' 


0. 9 


p 


2. 145966026 
1453. 383883 
2063 a 7516*28 
9422966. 381 


k 

k 

ka*' 

y ' 

Area 


. !;;i b 6 6 4 •' 1 6 y 


P 

k 


1 *..:' 5 4 . d 6 4 6 1 

1 QOC 0 7 77 • J £ 
1 7 u- 1 i 1 c. 

8184684. 605 


k o „ 
ko?; 

y 

Area 
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III. Program Listing 



Before entering the program, press 2nd and then CP or 
turn the calculator off and then on. Next enter 5 in the display, 
press 2nd and then Op 17. This repartitions the calculator's 
memory so that the complete program can be entered. 

Before recording the program, enter 6 in the display, 
press 2nd and then Op 17. This returns the calculator's memory 
to the normal partition (479.59). Returning the calculator to 
the normal partition allows the two program cards to be read 
in the normal partition without forcing. When the program 
is used, it repartitions the calculator so that Bank 3 
registers are program registers. 
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V 



000 


-—I •- 


LBL 


050 


47 


CMS 


001 


15 


E 


051 


15 


E 


no 2 


n 9 


9 


052 






003 


42 


STO 


053 


42 


STD 


no 4 


00 


00 


054 


09 


09 


005 


42 


STD 


055 


69 


DP 


fi 106 


01 


01 


056 


i”i 

i_i 


v Q 


007 


4 ? 


RTH 


057 


86 


STF 


008 


7 6 


LBL 


058 


00 


0 0 


009 


18 


C ' 


059 


32 


X i T 


0 1 0 


69 


DP 


060 


99 


PRT 


VJj 1 1 
_£LL2_ 
0 1 3 


20 

00 


20 o- 
ST* ^ 

00 


06 1 
062 
063 


18 

91 

99 


C 1 

R/S 

PRT 


014 


92 


RTH 


064 


is 


C ’ 


015 


,* t« 


LBL 


065 


91 


R/S 


0 1 6 


19 


n* 


066 


99 


PRT 


017 


69 


DP 


067 


18 


C ■ 


0 1 8 


21 


21 


063 


91 


R/S 


0 1 9 


I* O 


RC* — ' 


069 


99 


PRT 


020 


01 


01 


070 


in 


E 1 



021 


o_2 


RTH 


071 


18 


c ' 


022 


76 


LBL 


072 


O "7 

i 


IFF 


023 


10 


E * 


073 


03 


03 


024 


65 


x 


074 


01 


01 


‘025 


O Q 
i_i y 




075 


67 


67 


026 


cr cr 




076 


19 


n * 


027 


01 


1 


077 


r. -7 

i_* i 


IFF 


028 


08 


l~l 

1 


078 


01 


01 


029 


uu 


0 


07 y 


ou 


IJ IJ 


030 


95 


= 


080 


91 


91 


031 


92 


RTH 


08 1 


75 


- 


032 


7 6 


LBL 


082 


19 


H‘ 


033 


12 


B 


M o -Z* 


9F, 


= 


034 


O 

o o 


STF 


084 


94 


+ / - 


035 


01 


01 


035 


3 8 


S I H 


036 


76 


LBL 


086 


65 




037 


1.1 


ft 


087 


19 


D * 


038 


■”* 


X 1 T 


088 


61 


GTD 


039 


.-« 
O i 


IFF 


089 


01 


01 


040 


00 


00 


090 


02 


02 


041 


00 


00 


09 1 


85 


+ 


042 


59 


59 


092 


19 


D ' 


043 


O “7 

i 


IFF 


093 


32 


X 1 T 


044 


03 


03 


094 


19 


D 1 


045 


00 


00 


095 


•j - ”i 


INV 


046 


59 


59 


096 


•" -7 

o r 


P/R 


047 


05 


5 


097 


24 


CE 


048 


69 


DP 
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IV. A Development for the Procedure 

In the development for the estimation procedure given here, 
all angles are in radians and the assumptions stated in Section I 
apply. 

Figure 4 shows three bearing lines from the ith of n stations. 
One is the observed bearing line of an object. One of length fi ^ goe 

to the origin of an xy-coordinate 
system located at the object's un- 
known position. And one of length 
r^ goes to an initial estimate 
with known position but unknown 
coordinates (x,y) . Note, esti- 
mates for -x and -y estimate 
the object's position. To find 
estimates -x and -y, consider 

the arc coordinates u. = fi . (0. —<}>.) 

1 i i l 

of the observed bearing line and 
v^ = of the bearing line 

to the point (x,y) . They are de- 
fined by the three bearing lines 
and the circle of radius fi ^ which 
goes through the object's position 
and which is centered on the station 
as shown in Figure 4 . 




FIGURE 4. Problem Geometry. 
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By defining (0^-0^) (all angles in radians), 

u . = v. + w.. Note, 0.-8. is known, but 8.-<J>. is not. 
111 'll it 

However, can be expressed in terms of x and y, and, to first 

order, = x cos 0^ - y sin 0.^; so, if tan (8^^) - ( 8^ — 4>^ ) for 
i = l,2,...,n, that is, if (x,y) is relatively near the object's 
position, u^ - + x cos 0^ - y sin 0^ for i = l,2,...,n. 

In this development, the initial estimate is assumed to be 
relatively close enough to the object's position so that the above 
approximation for u^ can be used and so that can be taken 

equal to r^. With this assumption, u^ is a function of: 0^, 

the observed value of a random quantity; the known parameters r^ 
and 8^; and the unknown parameters x and y. 

If a distribution for the 0^ can be specified, then a 
distribution for the U. can be determined and standard estimates 

l 

A /V 

x and y for x and y can be considered. Maximum likelihood 

estimates are discussed in this section. Each 0^ is taken to 

2 

be a normal random variable with mean 6. and variance e.. And 

T i l 

the n random variables 0^, i = 1,2,. . . ,n (one for each station) 
are taken to be independent. 

The likelihood for a sample ’ ' e n ; '' S t ^ ien 




and the likelihood for a corresponding sample u^ , , u^ is 
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where a. = r.e. (with e. in radians) since u. = r . ( 0 . -<b . ) . 

111 1 111 T i 

By definition, the maximum likelihood estimates of x 

A A 

and y are the estimates x and y which make L (u^ , , . . . , u fi ) 

a maximum. In this case, making L (u, , u~ , . . . , u ) a maximum is 

1 2 n 

2 2 A 
equivalent to making (u^/a^) a minimum. So, to find x 

A 

and y, solve the following two equations for x and y: 



3 ( £n L) 
3x 



0 and 



3Un L) 
3y 



The solutions are x = x and y = y, and x and y are the 
maximum likelihood estimates. With w. = r. (0.-3.) and 

i ill 

the conditions assumed above these two equations are linear 
equations in x and y. And, 



n 



I [w. + x cos 3 . - y sin 6.] (cos 3 • ) /a . 

j l l J l l i 



2 = 0 



and 



n 



[w^ + x cos 3^ - y sin 3^] (sin 3^)/^^ = 0 



And, in terms of the following quantities: 

2 2 2 
A = E(cos 3 • ) /o . , B = E(sin 3. cos 3 • ) /a . , 

l ' i l i ' i 

C = E (sin 2 3 .)/a 2 , D = E(w. cos 3.)/a? , 

1 1 1 11 

E = E (w. sin 3 • ) /a 2 , 

l l l 



27 



the equations are : 

Ax - By = -D 

Bx - Cy = -E . 

So the solutions are: 

x = (BE - CD)/ (AC - B 2 ) 

y = (AE - BD) / (AC - B 2 ) . 

A confidence region can be constructed about an estimated 
position. In order to indicate how this is done, a probability 
region about the true position will be considered first. 

A A 

Note, x and y are values of random variables. If a 
new set of bearings 0-^,0 2 ,...,0 n is observed (for a fixed initial 
estimate and object), in general, a new pair of values x and 

will be obtained. 

/\ /\ 

If X and Y represent these random variables, then 

A -in p 

X = l (W./a.) (B sin 6- - C cos 6 • ) 

(AC-B ) 1 1 1 1 1 



Y = i- 5 - l (W./a 2 ) (A sin & - B cos 3.) 

(AC-B'*) 1 1 1 1 1 

with W^ = r^ (0^-6^) . (W^ is the random distance intercepted 

along the ith arc between the bearing lines defined by 0^ and 

V) 
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Note, X and Y have a bivariate normal distribution, 
since they are a linear combination of the n normal random 
variables w i , W 2 ' * * * • w n ' or e< 3 u i va l ent ly °f the n normal 
random variables 2 ' * * * '® n * A l so E(W^) = r^(cf)^-g^). 

If 6^ = <J>^ for i = l,2,...,n, that is, if the initial 
estimate of the object's position is at the object's position, 

E(W^) = 0 for i = 1,2, ...,n. And, therefore, E(X) = 0 and 
E(Y) =0. So, in this case, the "location" of the bivariate 

A A 

normal distribution of a point (x, Y) , the random 
coordinates of the object's estimated position, is the same as 
that for the point (-X, -Y) and both are centered on the 

object's position. However, the "location" of the distribution 

/\ s\ 

of (-X,-Y) is independent of the location of the initial 

A A 

estimate when the coordinates (-X,-Y) refer to a coordinate 
system with origin at the initial estimate. This fact simplifies 
the establishment of a confidence region about the location of 
an estimated position. 

A region of minimum area for a given probability of con- 
tainment of an estimated position can be determined. The region 
is bounded by an ellipse which is centered on the object's position 
and whose axes lie along the axes of an x ' y ' -coordinate system 
obtained by rotating the xy-coordinate system centered on the 

object's position through an angle y . In this system, = 0, 

~ * xy 

that is, X' and Y' are independent normal random variables. 
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The two coordinate systems are illustrated in Figure 5 . 
The coordinates of a point in the two systems are related by 



y ' y. 




FIGURE 5. Rotation Geometry. 



x' = x cos y + Y si n Y 
y' = -x sin Y + y cos y 

These relations, along with 

o ~ * =0, imply: 

x ' y ' 



2 2 2 2 2 
a = a cos y + 2a cos y sin Y + ° sin Y r 

A A A A A 

x ' x xy y 

2 2 

2 2 2 
a = a sin y “ 2a cos y sin y + a cos Y 

A A 1 A A 1 • A 1 

y ' x xy y 



tan 2 y = 2 X ^ ' 2 

o^-o 

A A 

x y 

where y, the angle of rotation of the coordinate axes, is 

positive in the counterclockwise direction. 

With the initial estimate of the object's position at 

the object's position (6 i = i = l,2,...,n), so E(W^) = 0 

and Var (W. ) = a? , 
l i ' 
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1 



o 2 = 



x (AC-B 2 ) 2 1 



r 2 2 

l (1/of) (B sin 6. - C cos 8.) , 



and 



o 2 = ~T~y I (l/o?) (A sin - B cos 8.) 2 

y (AC-B ) 1 1 1 1 

i n 2 

— T~ 7 1 (l/ a ,- ) (B sin 8 . - C cos 8 • ) (A sin 8 . - B cos 8 . ) 
-B ) 1 1 1 1 1 1 



0 /\ /\ *” 

xy (AC-B“) 

Using the definition for A, B and C, the above become 



J2 _ C 

a - 2 ' 

x (AC-B ) 

„2 _ A 

a ~ 2 ' 

y (AC-B ) 

and 

a = B 
xy (AC-B 2 ) 

So, tan 2y = 2B/(C-A) for 8^ = <t>^, i = l,2,...,n. 

With the object's position known and, hence, <Jk known 

2 2 

for i = l,2,...,n, the above equations for o , c , a and y 

A A A A 

x y xy 

can be used, since the initial estimate of the object's position 

can be taken as the object's position. 

With values for a A , a A , and y, values for 

x y xy x ' 

and can be found by using the equations in the middle of 

y' 

Page 30. And then, the probability that an estimated 
position will be within an ellipse of semiaxes ka^ and ka^ 
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which is centered on the object's position can be found. It is 
2 

1 - exp(-k /2) . (This result follows from integrating the bi- 
variate normal density over the ellipse.) And the area of the 

ellipse is 7 rk 2 a^ a * 

x ' y ' 

A, A 

Given estimates x and y found by using the relations 

on Page 28, an ellipse with semi-axes kg A and ka^ centered 

x' y ' 

on the point with coordinates (-x,-y) in a coordinate system 

with origin at the initial estimate and oriented as indicated 

2 

by y is a 1 - exp(-k /2) confidence region. This follows 

/N A 

from the bivariate normal distribution of -X and -Y which 

in this system is centered on the object's position. The ellipse 

2 2 

is defined if and a are known (the covariance matrix 

x y xy 

is known) . And to the degree of the approximations involved, 

this can be assumed to be the case.. In particular, by assuming 

the initial estimate of the object's position is at the object's 

position, which is consistent with assuming (6^-<J>^) is small, 

2 2 

values for o A , a A and Y can be obtained by using 

x y xy 

the relations on Page 31. These values can then be used to 
2 2 

determine and o A by using the relations on Page 30. 

x' y' 

And, then, with a value for k, a confidence region can be 
constructed. To the degree of the approximations involved, the 
shape of the confidence region is independent of both the object's 
position and of the initial estimate of the object's position. 
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For the case where bearings are taken from the object 
on two or more stations, e i is the reciprocal of the 
bearing taken from the object. 

A discussion for this and for other bearings only 
position estimation procedures for situations similar to the 
one considered here is given in Reference 1 listed below. 
Reference 2 gives an equivalent bearings only procedure. It 
also gives a range only procedure, a range and bearing procedure 
and HP-9830A programs with which to implement the procedures. 
Using the fix determined by two lines of bearing as the initial 
estimate was suggested by this reference. 



The equations used in the program to determine (x*, y*) , 
the coordinates of the fix, are: 



x* sin 
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